Edge effects in Bilayer Graphene Nanoribbons 
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We show that the ground state of zigzag bilayer graphene nanoribbons is non- magnetic. It also 
possesses a finite gap, which has a non-monotonic dependence with the width as a consequence of 
the competition between bulk and strongly attractive edge interactions. All results were obtained 
using ab initio total energy density functional theory calculations with the inclusion of parametrized 
van der Waals interactions. 
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Since the synthesis of graphene [1], a plethora of in- 
triguing properties has been found in this two dimen- 
sional zero gap crystal due to the presence of massless 
fermions with a high mobility [I, 0]. Besides the mono- 
layer, stacking two layers of graphene still preserves the 
high mobility, and some features of the electronic spec- 
trum can be controlled, for example, by applying an ex- 
ternal electric field [4]. Measurements of quantum hall 
effect and quasi-particle band structure indicate qualita- 
tive differences between monolayers and bi-layers. The 
occurrence of this rich physics at room-temperature [f| Q 
has attracted a great interest in designing graphene- 
based nanoelectronic devices. In this scenario, it is fun- 
damental to establish and control an energy gap (E g ). 

A possibility is to introduce lateral quantum confine- 
ment via synthesis of single layer (GNR)[7, 8] or bilayer 
(B-GNR)[9|, graphene nanoribbons by plasma etch- 
ing or chemical routes pj], EH- This opens a gap that 
rises the possibilities to use graphene in nanoelectronics, 
where small widths (sub- 10 nm) is required for room tem- 
perature applications [6]. However, the BGRNs are less 
sensitive to external perturbations in comparison with 
GNRs, and hence, they may be more appropriate to fab- 
ricate hight quality nano-devices [l3| . 

The electronic structure of the nanoribbons, including 
the gap, are largely affected by the geometric pattern 
(zigzag or armchair) at their edges. GNRs with zigzag 
edges, in particular, have as a distinct feature the pres- 
ence of edge states that introduce a large density of states 
(DOS) at the Fermi energy. Theoretical works predict 
that this configuration is unstable, and there will be the 
appearance of a magnetic order that leads to a removal 
of this large DOS peak@,[l,[uj. Magnetism in GNRs has 
been intensively investigated as a possible way to develop 
spintronic devices [HI, 0, Q3] • An anti- ferromagnetic (fer- 
romagnetic) order between the two edges leads to a semi- 
conductor (metallic) state @, H, 03] • It is a l so believed 
that magnetism is necessary to open a gap in zigzag B- 
GNR (B-ZGNR)0. 

In this Letter we investigate the geometrical and elec- 
tronic structure of B-ZGNR, and show that the ground 
state of these systems is non-magnetic and possesses a 



non-monotonic finite gap. There are two possible edge 
alignments for the B-ZGNR, called a and (3 (see Fig. 
1). We found that the a alignment is energetically fa- 
vorable, with an inter-layer edges attraction, whereas 
for the /3 alignment there is an inter-layer edges repul- 
sion. These edge-related forces cause a deviation from 
the exact Bernal stacking, resulting in a non-monotonic 
behavior of the energy gap with the width w for the a 
B-ZGNR, with a maximum value at w ~ 3.5nm. These 
results differ qualitatively from GNRs with zigzag edges 
(M-ZGNR)[ll|,[l3 

All our results are based on ab initio total energy Den- 
sity Functional Theory [19] (DFT) calculations. In order 
to correctly describe multi-layer graphitic compounds, it 
is necessary to include van der Waals (vdW) interactions. 
The use of fully relaxed total energy DFT calculations to 
study such systems suffers from serious limitations, since 
the most traditional exchange correlation (xc) function- 
al in use today do not correctly describe these terms. 
With the LDA xc, the geometry is correctly described 
but the inter-layer binding energy is underestimated by 
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FIG. 1: Bilayer graphene nanoribbons with (a) a and (b) (3 
edge alignment, (c) Side view of bilayer graphene nanorib- 
bons. The blue (white) atoms form the upper (bottom) layer. 



TABLE I: Inter-layer binding energy (Eb), in eV/atom, and 
distance h, in A, for graphite and a graphene bilayer, which 
is not bound (NB) at the GGA level. 
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50%, whereas the GGAs xc do not even correctly de- 
scribe the geometrical features [20]. Thus, in order to be 
able to investigate the geometries and relative energies of 
B-ZGNRs, we include a non-local potential in the Kohn- 
Sham (KS) equations that correctly describes the vdW 
interactions. We modified the SIESTA code[23[ adding 
in the KS Hamiltonian[24] the dispersion corrected atom 
centered potential (DCACP)[25]. This correction is suf- 
ficiently accurate to describe weakly bonded systems (26| 
with the vdW interactions included in the whole self- 
consist cycle, providing accurate values for both forces 
(and thus geometries) and total energies. Our imple- 
mentation was successfully tested (see Table |B, and was 
employed to obtain the results here reported [28i l29l. l30j . 

We investigated B-ZGNR composed by two M-ZGNR 
passivated with hydrogen, and widths [31] that range 
from w — 0.6 to w = 4.5 nrn. The layers are in the 
Bernal stacking, which means that there are two types of 
C atoms, those that are positioned above the center of 
the hexagons of the other layer, defining a B-sublattice, 
and those right on top of the C atoms of the other layer, 
forming an A-sublattice. An infinite graphene bilayer has 
no gap, and the orbit als at the Fermi level are located at 
the B-sublattice. When we cut the layer along the zigzag 
edge, there are two possible alignments (FigHJ): (a) the a 
alignment, where the outermost edge atoms belong to the 
A-sublattice, and (b) the (3 alignment, where the outer- 
most edge atoms belong to the B-sublattice. Thus, only 
the inter-layer edge interaction differs. Two geometri- 
cal distortions have proven to be important: (i) an edge 
distortion that causes a curvature in the ribbons (see 
Fig life)), and (ii) a lateral deviation from the perfect 
Bernal stacking. To quantify this deviation, we define 
the quantity u = dc-c — d, where d is shown in FigQJa) 
and (b), and dc-c is the carbon-carbon bond length. 
The perfect Bernal stacking corresponds to u = 0. 

The geometries and band structures of fully relaxed 
B-ZGNRs with w ~ 1.0 nm for a and (3 alignments are 
presented in Fig. O In B-ZGNR with the (3 alignment, 
similarly to M-ZGNRs, a non spin-polarized calculations 
leads to a high DOS at the Fermi energy, and a magnetic 
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FIG. 2: Ground State of fully relaxed B-ZGNRs generated by 
stacking two (5,0) M-ZGNR. Below each band structure the 
geometry and local magnetization are presented, (a) a align- 
ment. This state is non-magnetic and presents a geometric 
distortion near the edge, (b) (3 alignment. This state shows 
an AF in-layer and AF inter-layer magnetic order. 



order is required to split these localized edge states at 
the K symmetry point. In order to establish the pos- 
sible spin polarized configurations, we used four initial 
guesses for the density matrix before starting the self- 
consistency cycle, which are: i) anti-ferromagnetic (AF) 
in-layer and inter-layer, ii) ferromagnetic (F) in-layer and 
inter- layer, Hi) AF in- layer and F inter- layer and iv) F in- 
layer and AF inter-layer. As well as non-polarized calcu- 
lations. From all calculations, the AF in-layer and inter- 
layer guess leads to the lower energy state (Fig. Hfb)). 
However, the energy differences are less than /crT[18|. 

At the a alignment (Fig. [2(a)), on the other hand, we 
obtain a qualitatively new situation. There is a strong 
attractive interaction between the edge atoms of the two 
layers, with a resulting geometric distortion that de- 
creases the distance between them (Fig. [2(a)). For all a 
B-ZGNRs that we have investigated, the final geometry 
always had an inter-layer edge atoms distance around 3.0 
A. The final configuration is non-magnetic and with a fi- 
nite gap, contrary to previous results where the presence 
of a gap was intrinsically coupled to a magnetic state [Tit. 
Note that if we do not allow the atoms at the two layers to 
relax, but simply optimize the inter-layer distance (i.e., 
the layers keep their planar geometries), a magnetic con- 
figuration is still necessary to open a gap[l8|. However, 
this configuration has higher energy. 

If we take one of the mono-layers that form the final 
relaxed a B-ZGNR, and perform a calculation without 
letting the atoms relax, we obtain an energy increase, 
when compared to the lowest energy M-ZGNR (AF in- 
layer [15]), that can be broken down in two components. 
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Considering M-ZGNR with widths larger than 2 nm, if 
we allow the distorted monolayer to be magnetic, the en- 
ergy increase is ~ 0.1 eV/nm, which can be viewed as the 
elastic contribution. If we now consider a non-magnetic 
configuration for this distorted monolayer, which is the 
situation in the a B-ZGNR, there is an extra energy in- 
crease of ~ 0.4 eV/nm, i.e., an overall energy penalty of 
^0.5 eV/nm. Considering the two monolayers, the total 
energy cost to deform and demagnetize the a B-ZGNR is 
~ 0.96 eV/nm. This energy increase is more than com- 
pensated by the edge atoms interaction, and the energy 
gain is in part associated with a large split of the localized 
edge states. At the f3 B-ZGNR (Fig. [2(b)), the unique 
way to diminish the DOS associated with the localized 
edge states at the Fermi energy is via a magnetic order- 
ing, and the system tends to increase the inter-layer edge 
atoms distance in order to allow a bigger magnetization, 
given rise to a repulsive inter-layer edge interaction. 
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FIG. 3: Dependence of the binding energies with the width. 

Comparing the two alignments, the a B-ZGNR results 
energetically favorable. This is an even more important 
conclusion conside ring that most of the calculations used 
the B-ZGNR[To|,[32|. Fig. [3]presents the dependence of 
the B-ZGNRs binding energies, for both edge alignments, 
as a function of w (calculated relative to two isolated low- 
est energy M-ZGNR). The interaction between the layers 
can be separated into two components: i) edge interac- 
tions, that do not depends on the width, ii) and bulk 
interactions, that increase linearly with w. The binding 
energies (per unit length) can be well adjusted with: 



Ejj(w) = a + bw. 



(i) 



Since there are two edges, a/2 is the inter-layer edge in- 
teraction energy per unit length, and b is the inter-layer 
bulk interaction energy per unit area. At the a align- 
ment, a = —0.26eV/nm and b = —2.0eV/nm 2 , indicating 
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FIG. 4: Dependence of the (a) energy gap (inset illustrates the 
change of character of the VBM) , and (b) the lateral deviation 
u with the width w (inset indicates how u affects the inter- 
layer A-sublattice p-orbitals interaction) . 



that the edges interaction is attractive (a < 0). For the (3 
alignment, a = +0.13eV/nm and b = —2.0eV/nm 2 , indi- 
cating that there is a repulsive edges interaction (a > 0), 
showing that its stability results solely from the bulk. 
The parameter 6, as expected, does not depend on the 
edge alignment, and it is very close to our calculated 
bulk inter-layer interaction in a graphene bilayer (0.027 
eVj atom=— 1 .99eV/ nm 2 ). 

For the a B-ZGNR, there is a competition between the 
forces deriving from the bulk, that do prefer the Bernal 
pattern of stacking, and the forces deriving from the 
edges, that tend to maintain the inter-layer edge carbon 
atoms distance close to 3 A. The system then minimizes 
the overall energy penalty by simultaneously optimizing 
both the deviation u from the exact Bernal stacking, and 
the elastic energy associated with the ribbon's curvature. 
Thus, for narrower B-ZGNR the system prefers to have a 
larger value of u and a smaller overall curvature whereas 
for wider B-ZGNR the deviation u tends to decrease to 
minimize the bulk energy penalty, since now it is possible 
to have a softer curvature that is somewhat localized at 
the edges when compared to the total width of the rib- 
bon (see Fig. [2(c)). As a result, the dependence of the 
lateral deviation u with the width w is well adjusted by 
u(w) = 0.80e~ - 35w (with u in [A], and w in [nm]). More- 
over, as a result, the carbon-carbon bond lengths do not 
significantly differ from their values in the M-ZGNRs. 
For this situation, the average inter-layer distance h is 
close to its value in the graphene bilayer (see Tab. [Ij). 

Without the geometrical deformation caused by the 
inter-layer edge interactions, a monotonic decrease of the 
energy gap is expected due to the quantum confinement 
(oc 1/w) [8l. Il8j. However, for small ribbons, we find that 
the character of the valence band maximum (VBM) is 
located at the A-sublattice, as opposed to larger rib- 
bons where it is located in the B-sublattice (see Fig. 
Ufa)), similarly to the Fermi level orbit als in the infi- 
nite graphene bilayer. Moreover, since the interaction 
between the C atoms in the A-sublattice increases when 
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u decreases (see FigHJb)), the gap initially increases with 
w. However, for w > 3.5 nm, due to the quantum- 
confinement decrease, there is a crossover between the 
two highest occupied bands, and the character of the 
VBM is at the B-sublattice. Thus, this leads to a non- 
monotonic behavior of the energy gap with w (Fig. Ufa)). 

For the /3 B-ZGNR, there is a repulsive interaction 
between the edges, in such way that the inter-layer edge 
carbon atoms distance is close to 3.7 A. There happens 
a small negative lateral deviation (u < 0) that can be 
neglected when w > 1.6nm. Note that for the a B- 
ZGNR, due to the attractive edges interaction, u > 0. 
We also found that, despite the presence of a magnetic 
order, the energy gap disappears when w > 3 nm. 

Summarizing, we unequivocally show that for B-ZGNR 
the edge alignment a is the lowest energy configuration. 
This is a result of the strong attractive interaction be- 
tween the edges, that is manifested in an observed chem- 
ical bonding between the inter-layer edge carbon atoms, 
and that significantly influences the geometry and elec- 
tronic structure of bilayer nanoribbons with sub- 10 nm 
widths. As a consequence, the ground- state is non- 
magnetic and possesses a finite gap, which presents a 
non-monotonic dependence with the width. 
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